4. Linear Response 

The goal of response theory is to figure out how a system reacts to outside influences. 
These outside influences are things like applied electric and magnetic fields, or applied 
pressure, or an applied driving force due to some guy sticking a spoon into a quantum 
liquid and stirring. 

We've already looked at a number of situations like this earlier in the course. If 
you apply a shearing force to a fluid, its response is to move; how much it moves is 
determined by the viscosity. If you apply a temperature gradient, the response is for 
heat to flow; the amount of heat is determined by the thermal conductivity. However, in 
both of these cases, the outside influence was time independent. Our purpose here is to 
explore the more general case of time dependent influences. As we'll see, by studying 
the response of the system at different frequencies, we learn important information 
about what's going on inside the system itself. 

4.1 Response Functions 

Until now, our discussion in this course has been almost entirely classical. Here we 
want to deal with both classical and quantum worlds. For both cases, we start by 
explaining mathematically what is meant by an outside influence on a system. 

Forces in Classical Dynamics 

Consider a simple dynamical system with some generalized coordinates Xi(t) which 
depend on time. If left alone, these coordinates will obey some equations of motion, 

'±i + gi{x,x) = 

This dynamics need not necessarily be Hamiltonian. Indeed, often we'll be interested 
in situations with friction. The outside influence in this example arises from perturbing 
the system by the addition of some driving forces Fi(t), so that the equations of motion 
become, 

x i + g i (x,x) = F i (t) (4.1) 

In this expression, Xi(t) are dynamical degrees of freedom. This is what we're solving 
for. In contrast, Fi(t) are not dynamical: they're forces that are under our control, like 
someone pulling on the end of a spring. We get to decide on the time dependence of 
each Fiit). 
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It may be useful to have an even more concrete example at the back of our minds. 
For this, we take every physicist's favorite toy: the simple harmonic oscillator. Here 
we'll include a friction term, proportional to 7, so that we have the damped harmonic 
oscillator with equation of motion 



We will discuss this model in some detail in section 4.2. 
Sources in Quantum Mechanics 

In quantum mechanics, we introduce the outside influences in a slightly different man- 
ner. The observables of the system are now operators, 0,. We'll work in the Heisenberg 
picture, so that the operators are time dependent: O = 0(t). Left alone, the dynamics 
of these operators will be governed by a Hamiltonian H(0). However, we have no 
interest in leaving the system alone. We want to give it a kick. Mathematically this is 
achieved by adding an extra term to the Hamiltonian, 



The (f>i(x) are referred to as sources. They are external fields that are under our 
control, analogous to the driving forces in the example above. Indeed, if we take a 
classical Hamiltonian and add a term of the form x<f> then the resulting Euler-Lagrange 
equations include the source on the right-hand-side in the same way that the force 
F appears in (4.2). 

4.1.1 Linear Response 

We want to understand how our system reacts to the presence of the source or the 
driving force. To be concrete, we'll chose to work in the language of quantum mechanics, 
but everything that we discuss in this section will also carry over to classical systems. 
Our goal is to understand how the correlation functions of the theory change when we 
turn on a source (or sources) </>i(x). 

In general, it's a difficult question to understand how the theory is deformed by the 
sources. To figure this out, we really just need to sit down and solve the theory all over 
again. However, we can make progress under the assumption that the source is a small 
perturbation of the original system. This is fairly restrictive but it's the simplest place 
where we can make progress so, from now on, we focus on this limit. Mathematically, 
this means that we assume that the change in the expectation value of any operator is 
linear in the perturbing source. We write 



x + 7± + UqX = F(t) 



(4.2) 



H sourcc (t) = 4>i(t)Oi(t) 



(4.3) 




(4.4) 
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Here Xij(t\ t') is known as a response function. We could write a similar expression for 
the classical dynamical system (4.1), where S(Oi) is replaced by Xi(t) and <fi is replaced 
by the driving force Fj(t). In classical mechanics, it is clear from the form of the 
equation of motion (4.1) that the response function is simply the Green's function for 
the system. For this reason, the response functions are often called Green's functions 
and you'll often see them denoted as G instead of x- 

From now on, we'll assume that our system is invariant under time translations. In 
this case, we have 

Xij(t;lf)=Xij(t-lf) 

and it is useful to perform a Fourier transform to work in frequency space. We define 
the Fourier transform of the function f(t) to be 

f(u) = J dte^f(t) and f(t) = J ^V^/M (4.5) 

In particular, we will use the convention where the two functions are distinguished only 
by their argument. 

Taking the Fourier transform of (4.4) gives 

6{Oi(u))= I dt' fdte^Xijit-t)^) 



= j dt' j dte^-^xdt-t') e^'0#) 

= xa(u) M u ) ( 4 - 6 ) 

We learn the response is "local" in frequency space: if you shake something at frequency 
iv, it responds at frequency iv. Anything beyond this lies within the domain of non- 
linear response. 

In this section we'll describe some of the properties of the response function xi 00 ) 
and how to interpret them. Many of these properties follow from very simple physical 
input. To avoid clutter, we'll mostly drop both the i,j indices. When there's something 
interesting to say, we'll put them back in. 

4.1.2 Analyticity and Causality 

If we work with a real source and a Hermitian operator O (which means a real 
expectation value (O)) then x(t) must also be real. Let's see what this means for the 
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Fourier transform x( u )- useful to introduce some new notation for the real and 
imaginary parts, 



This notation in terms of primes is fairly odd the first time you see it, but it's standard 
in the literature. You just have to remember that, in this context, primes do not mean 
derivatives! 

The real and imaginary parts of the response function x{°°) have different interpre- 
tations. Let's look at these in turn 

• Imaginary Part: We can write the imaginary piece as 



We see that the imaginary part of x( u ) * s due to the part of the response func- 
tion that is not invariant under time reversal t — > —t. In other words, x"^) 
knows about the arrow of time. Since microscopic systems are typically invariant 
under time reversal, the imaginary part x"^) must be arising due to dissipative 
processes. 

X"(w) is called the dissipative or absorptive part of the response function. It is 
also known as the spectral function. It will turn out to contain information about 
the density of states in the system that take part in absorptive processes. We'll 
see this more clearly in an example shortly. 

Finally, notice that x"^) is an °dd function, 



x>) = ~M")-x»] 






• Real Part: The same analysis as above shows that 
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The real part doesn't care about the arrow of time. It is called the reactive part 
of the response function. It is an even function, 



x'(-oj) = +x'(oj) 



Before we move on, we need to briefly mention what happens when we put the labels 
i, j back on the response functions. In this similar analysis to that above shows 

that the dissipative response function comes from the anti-Hermitian part, 



Causality 

We can't affect the past. This statement of causality means that any response function 
must satisfy 



For this reason, x is often referred to as the causal Green's function or retarded Green's 
function and is sometimes denoted as Gr(£). Let's see what this simple causality 
requirement means for the Fourier expansion of Xi 



When t < 0, we can perform the integral by completing the contour in the upper-half 
place (so that the exponent becomes — iu x (— i\t\) — > — oo). The answer has to be zero. 

Of course, the integral is given by the sum of the residues inside the contour. So if 
we want the response function to vanish for all t < 0, it must be that x( u ) h &s no poles 
in the upper-half plane. In other words, causality requires: 



4.1.3 Kramers-Kronig Relation 

The fact that x is analytic in the upper-half plane means that there is a relationship 
between the real and imaginary parts, x' an d x" ■ This is called the Kramers-Kronig 
relation. Our task in this section is to derive it. We start by providing a few general 
mathematical statements about complex integrals. 



1 



(4.7) 



X(t) = for all t < 




x(w) is analytic for Imu > 
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A Discontinuous Function 

First, consider a general function p(ou). We'll ask that p(u) is meromorphic, meaning 
that it is analytic apart from at isolated poles. But, for now, we won't place any 
restrictions on the position of these poles. (We will shortly replace p(uo) by x(°°) which, 
as we've just seen, has no poles in the upper half plane). We can define a new function 
f{uj) by the integral, 

f M = lf b J^l^ (4.8) 



Here the integral is taken along the interval u/ G [a, b] of the real line. However, when 
uj also lies in this interval, we have a problem because the integral diverges at uj' = uj. 
To avoid this, we can simply deform the contour of the integral into the complex plane, 
either running just above the singularity along uj' + ie or just below the singularity 
along uj' — ie. Alternatively (in fact, equivalently) we could just shift the position of 
the singularity to uj — > uj =)= e. In both cases we just skim by the singularity and the 
integral is well defined. The only problem is that we get different answers depending 
on which way we do things. Indeed, the difference between the two answers is given by 
Cauchy's residue theorem, 

l -[f(oo + ie)-f(uo-ie)]=p(uj) (4.9) 

The difference between f(uj+ie) and f(uj—ie) means that the function f(u) is discontin- 
uous across the real axis for uj e [a, b]. If p{yj) is everywhere analytic, this discontinuity 
is a branch cut. 

We can also define the average of the two functions either side of the discontinuity. 
This is usually called the principal value, and is denoted by adding the symbol V before 
the integral, 

\[f(uj + ie) + f(uj - ie)} = -V f 4^-duj' (4.10) 
2 in J a uj — uj 

We can get a better handle on the meaning of this principal part if we look at the real 
and imaginary pieces of the denominator in the integrand 1/[uj' — (uj ± ie)], 

1 uj' — uj ie , 

+ ,<> , (4-11) 



uj'-(uj± ie) (uj' - uj) 2 + e 2 (uo' - uj) 2 + e 2 

By taking the sum of /(a; -Me) and f(uj — ie) in (4.10), we isolate the real part, the first 
term in (4.11). This is shown in the left-hand figure. It can be thought of as a suitably 
cut-off version of 1/(uj' — uj). It's as if we have deleted an small segment of this function 
lying symmetrically about divergent point uj and replaced it with a smooth function 
going through zero. This is the usual definition of the principal part of an integral. 
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Figure 9: The real part of the function Figure 10: The imaginary part of the 
(4.11), plotted with u' = 1 and e = 0.5. function (4.11), plotted with u/ = 1 and 

e = 0.5 



We can also see the meaning of the imaginary part of 1/ (u/ — to) , the second term 
in (4.11). This is shown in the right-hand figure. As e — > 0, it tends towards a delta 
function, as expected from (4.9). For finite e, it is a regularized version of the delta 
function. 

Kramers-Kronig 

Let's now apply this discussion to our response function xi^)- We'll be interested in 
the integral 



where the contour C skims just above the real axis, before closing at infinity in the 
upper-half plane. We'll need to make one additional assumption: that x( z ) f& ns on? 
faster than l/\z\ at infinity If this holds, the integral is the same as we consider in 
(4.8) with [a, b] — > [— oo, +oo]. Indeed, in the language of the previous discussion, the 
integral is f{u — ie), with p = X- 

We apply the formulae (4.9) and (4.10). It gives 



But we know the integral in (4.12) has to be zero since x{ u ) nas no poles in the 
upper-half plane. This means that f(uj — ie) = 0, or 




(4.12) 





(4.13) 
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The important part for us is that factor of sitting in the denominator. Taking real 
and imaginary parts, we learn that 

Re X (co)=V 4.14 

and 

T / x ^ f + °° duj'Re X (uj') 

Imx(w) = -P / 4.15 

These are the Kramers- Kronig relations. They follow from causality alone and tell 
us that the dissipative, imaginary part of the response function x"^) is determined 
in terms of the reactive, real part, x'i 00 ) an d vice-versa. However, the relationship is 
not local in frequency space: you need to know x'i 00 ) f° r an frequencies in order to 
reconstruct x" f° r an y single frequency. 

There's another way of writing these relations which is also useful and tells us how 
we can reconstruct the full response function xi 00 ) if we onr y know the dissipative part. 
To see this, look at 

dJ Imx(a/) (4.16) 



m uj' — uj — ie 



where the — ie in the denominator tells us that this is an integral just below the real 
axis. Again using the formulae (4.9) and (4.10), we have 



Im X (u/) f + °°duj' Im X («/) 

Imx(w) + V 



oo m uj' — uj — ie 



= ImxH-iItexH (4.17) 
Or, rewriting as x( u ) = R- e x( w ) + 2lmx(u;), we get 



If you know the dissipative part of the response function, you know everything. 
An Application: Susceptibility 

Suppose that turning on a perturbation induces a response (O) for some observable 
of our system. Then the susceptibility is defined as 

d(0) 



X = 



u)=0 



- 85 - 



We've called the susceptibility x which is the same name that we gave to the response 
function. And, indeed, from the definition of linear response (4.4), the former is simply 
the zero frequency limit of the latter: 



A common example, which we met in our first course in statistical mechanics, is the 
change of magnetization M of a system in response to an external magnetic field B. 
The aptly named magnetic susceptibility is given by x — OM/dB. 

From (4.18), we can write the susceptibility as 



We see that if you can do an experiment to determine how much the system absorbs 
at all frequencies, then from this information you can determine the response of the 
system at zero frequency. This is known as the thermodynamic sum rule. 

4.2 Classical Examples 

The definitions and manipulations of the previous section can appear somewhat ab- 
stract the first time you encounter them. Some simple examples should shed some 
light. The main example we'll focus on is the same one that accompanies us through 
most of physics: the classical harmonic oscillator. 

4.2.1 The Damped Harmonic Oscillator 

The equation of motion governing the damped harmonic oscillator in the presence of a 
driving force is 



Here 7 is the friction. We denote the undamped frequency as ou , saving u for the 
frequency of the driving force as in the previous section.. We want to determine the 
response function, or Green's function, x(t — t') of this system. This is the function 
which effectively solves the dynamics for us, meaning that if someone tells us the driving 
force F(t), the motion is given by 




(4.19) 



x + 7X + UqX = F(t) 



(4.20) 




(4.21) 
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Figure 11: The real, reactive part of the response function for the under damped harmonic 
oscillator, plotted with loq = 1 and 7 = 0.5. 



There is a standard method to figure out Firstly, we introduce the (inverse) 

Fourier transform 

X(t) = / ^-xH 

We plug this into the equation of motion (4.20) to get 

"°du 



2tt 



dt'[-uo 2 - iju + ujfa-^-^xiuW) = F(t) 



which is solved if the J duo gives a delta function. But since we can write a delta 
function as 2n5(t) = J due~ luJt , that can be achieved by simply taking 

XM = ^ 2 ( 4 - 22 ) 

There's a whole lot of simple physics sitting in this equation which we'll now take some 
time to extract. All the lessons that we'll learn carry over to more complicated systems. 

Firstly, we can look at the susceptibility, meaning x{w = 0) = 1/u)q. This tells us 
how much the observable changes by a perturbation of the system, i.e. a static force: 
x = F/uq as expected. 

Let's look at the structure of the response function on the complex cu-plane. The 
poles sit at ul + i^u^ — uj^ = or, solving the quadratic, at 



There are two different regimes that we should consider separately, 
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Figure 12: The imaginary, dissipative part of the response function for the underdamped 
harmonic oscillator, plotted with luq = 1 and 7 = 0.5. 

• Underdamped: Uq > 7 2 /4. In this case, the poles have both a real and imag- 
inary part. They both sit on the lower half plane. This is in agreement with 
our general lesson of causality which tells us that the response function must be 
analytic in the upper-half plane 

• Overdamped: Uq < 7 2 /4. Now the poles lie on the negative imaginary axis. 
Again, there are none in the upper-half place, consistent with causality. 

We can gain some intuition by plotting the real and imaginary part of the response 
function for oj e R. Firstly, the real part is shown in Figure 11 where we plot 



This is the reactive part. The higher the function, the more the system will respond 
to a given frequency. Notice that Rexioj) is an even function, as expected. 

More interesting is the dissipative part of the response function, 



This is an odd function. In the underdamped case, this is plotted in Figure 12. Notice 
that Im% is proportional to 7, the coefficient of friction. The function peaks around 
±co>o> at frequencies where the system naturally vibrates. This is because this is where 
the system is able to absorb energy. However, as 7 — > 0, the imaginary part doesn't 
become zero: instead it tends towards two delta functions situated at ±uj . 




(4.23) 




(4.24) 
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4.2.2 Dissipation 

We can see directly how lmx{uj) is related to dissipation by computing the energy 
absorbed by the system. This what we used to call the work done on the system before 
we became all sophisticated and grown-up. It is 

— = F(tHt) 

+OC 



F(t)j t J_yt'x(t-t')F(t') 

/ + OO p+OO J 

dt' / ^{-i^e-^'hi^Fit') 

I 



—00 

p^i-iuxi^e-^^F^F^') (4.25) 

Let's drive the system with a force of a specific frequency Q, so that 

F(t) = F cos Vlt = F Re(e~ int ) 

Notice that it's crucial to make sure that the force is real at this stage of the calculation 
because the reality of the force (or source) was the starting point for our discussion 
of the analytic properties of response functions in section 4.1.2. In a more pedestrian 
fashion, we can see that it's going to be important because our equation above is not 
linear in F{u), so it's necessary to take the real part before progressing. Taking the 
Fourier transform, the driving force is 

F(u) = 2ttF [8(u + 5(uj + n)] 

Inserting this into (4.25) gives 

^ = -tFfa [ X (^)e- int - X (-^)e +tnt ] [e~ int + e m ] (4.26) 

This is still oscillating with time. It's more useful to take an average over a cycle, 
dW Q f 2n/n , dW 



dt 



= -jf dt— = -iFfa[ x (n)- x (-m 



But we've already seen that Rex(^) is an even function, while Imx(w) is an odd 
function. This allows us to write 



^ = 2F 2 fiIm X (fi) (4.27) 
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We see that the work done is proportional to Imx- To derive this result, we didn't 
need the exact form of the response function; only the even/odd property of the 
real/imaginary parts, which follow on general grounds. For our damped harmonic 
oscillator, we can now use the explicit form (4.24) to derive 



dW 2p2 7^ 



dt u (wg - Q 2 ) 2 + (7O) 2 

This is a maximum when we shake the harmonic oscillator at its natural frequency, 
Q = uo. As this example illustrates, the imaginary part of the response function tells 
us the frequencies at which the system naturally vibrates. These are the frequencies 
where the system can absorb energy when shaken. 



4.2.3 Hydrodynamic Response 

For our final classical example, we'll briefly return to the topic of hydrodynamics. One 
difference with our present discussion is that the dynamical variables are now functions 
of both space and time. A typical example that we'll focus on here is the mass density, 
p(x,t). Similarly, the driving force (or, in the context of quantum mechanics, the 
source) is similarly a function of space and time. 



Rather than playing at the full Navier-Stokes equation, here we'll instead just look 
at a simple model of diffusion. The continuity equation is 

We'll write down a simple model for the current, 

j=-DVp + F (4.28) 



where D is the diffusion constant and the first term gives rise to Fick's law that we met 
already in Section 1. The second term, F = F(x,t), is the driving force. . Combining 
this with the continuity equation gives, 

^ - DV 2 p — — V • F (4.29) 
We want to understand the response functions associated to this force. This includes 

— * 

both the response of p and the response of J, 
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For simplicity, let's work in a single spatial dimension so that we can drop the vector 
indices. We write 



where we've called the second label J on both of these functions to reflect the fact that 
F is a driving force for J. We follow our discussion of Section 4.1.1. We now assume 
that our system is invariant under both time and space translations which ensures that 
the response function depend only on t' — t and x' — x. We then Fourier transform with 
respect to both time and space. For example, 

*m = 

Then in momentum and frequency space, the response functions become 

p(u,k) = x P j(u,k)F(uj,k) 
J(w,k) = xjj(u,k)F(u,k) 

The diffusion equation (4.29) immediately gives an expression for XpJ- Substituting 
the resulting expression into (4.28) then gives us x.jj- The response functions ar 

ik —ioj 



-iu + Dk 2 ' /XJJ -iu + Dk 2 

Both of the denominator have poles on the imaginary axis at uo = —iDk 2 . This is the 
characteristic behaviour of response functions capturing diffusion. 

Our study of hydrodynamics in Sections 2.4 and 2.5 revealed a different method of 
transport, namely sound. For the ideal fluid of Section 2.4, the sound waves travelled 
without dissipation. The associated response function has the form 

1 

Xsound 



uj 2 — v 2 k 2 



which is simply the Green's function for the wave equation. If one includes the effect of 
dissipation, the poles of the response function pick up a (negative) imaginary part. For 
sound waves in the Navier-Stokes equation, we computed the location of these poles in 
(2.76). 
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4.3 Quantum Mechanics and the Kubo Formula 

Let's now return to quantum mechanics. Recall the basic set up: working in the 
Heisenberg picture, we add to a Hamiltonian the perturbation 

^sourcc(t) =0i(t)0i(t) (4.30) 

where there is an implicit sum over j, labelling the operators in the theory and, corre- 
spondingly, the different sources that we can turn on. Usually in any given situation we 
only turn on a source for a single operator, but we may be interested in how this source 
affects the expectation value of any other operator in the theory, {Oj). However, if we 
restrict to small values of the source, we can address this using standard perturbation 
theory. We introduce the time evolution operator, 

[/(Mo) = Texp H SOUIce (t')dt^j 

which is constructed to obey the operator equation idU jdt = H source U. Then, switching 
to the interaction picture, states evolve as 

m))i = u(t,t )mto))i 

We'll usually be working in an ensemble of states described by a density matrix p. If, 
in the distant past t — > oo, the density matrix is given by po, then at some finite time 
it evolves as 

pit) = U(t)p Q U-\t) 

with U(t) = U(t,t — > — oo). From this we can compute the expectation value of any 
operator Oj in the presence of the sources 0. Working to first order in perturbation 
theory (from the third line below), we have 

= Tr PQ (t)U-\t)Oi(t)U(t) 

« Tr p (t) (Oi{t) + i j dt' [H somce (t'), Oi(t)] + 
= (Oi{t)% =Q + i f dt'([H SOUTCC (t'),O t (t)}) + ... 

J — oo 

Inserting our explicit expression for the source Hamiltonian gives the change in the 
expectation value, 5(Oi) = (Oi)^ — (Cj)^=o, 

S(Oi) = i f dt' {[0 3 {t'),0 % {t)])^it') 

/+oo 
dt'e{t-t'){\0 3 {t'),O l {t)]) < P 3 {t') (4.31) 
-oo 
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where, in the second line, we have done nothing more than use the step function to 
extend the range of the time integration to +00. Comparing this to our initial definition 
given in (4.4), we see that the response function in a quantum theory is given by the 
two-pont function, 

Xij (t - t>) = -i6{t - t>) {[Oi(t), Oj(t')\) (4.32) 

This important result is known as the Kubo formula. (Although sometimes the name 
"Kubo formula" is restricted to specific examples of this equation which govern trans- 
port properties in quantum field theory. We will derive these examples in Section 4.4). 

4.3.1 Dissipation Again 

Before we make use of the Kubo formula, we will first return to the question of dis- 
sipation. Here we repeat the calculation of 4.2.2 where we showed that, for classical 
systems, the energy absorbed by a system is proportional to Im%. Here we do the same 
for quantum systems. The calculation is a little tedious, but worth ploughing through. 

As in the classical context, the work done is associated to the change in the energy 
of the system which, this time, can be written as 

dW d 

_ = _T r pH = T r (pH + pH) 

To compute physical observables, it doesn't matter if we work in the Heisenberg or 
Schrodinger picture. So lets revert momentarily back to the Schrodinger picture. Here, 
the density matrix evolves as ip = [H,p], so the first term above vanishes. Meanwhile, 
the Hamiltonian H changes because we're sitting there playing around with the source 
(4.30), providing an explicit time dependence. To simplify our life, we'll assume that 
we turn on just a single source, 0. Then, in the Schrodinger picture 

H = 0<p{t) 

This gives us the energy lost by the system, 
dW 

— = tt( p o 0) = (o)tf = [<e% =0 + 5(0)\ 

We again look at a periodically varying source which we write as 

0(t) = Re(0 o e- im ) 
and we again compute the average work done over a complete cycle 

HW VI f 2w/n , dW 
dt 



dt 2n J Q dt 
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The term (O(x))o cancels out when integrated over the full cycle. This leaves us with 

~dW Q f - : " , , , , , 

dt I dt' x(t ~ t ) <j>(t ) <j>(t) 



dt 



2tt 



oo 

q i>2ir/Q p+oo i 



x 



where the 2 and 0* 2 terms have canceled out after performing the J dt. Continuing, 
we only need the fact that the real part of x is even while the imaginary part is odd. 
This gives us the result 



dW 
~df 



2fi X "(fi)|0 o | 



(4.33) 



Finally, this calculation tells us about another property of the response function. If we 
perform work on a system, the energy should increase. This translates into a positivity 
requirement > 0. More generally, the requirement is that Vtx'ljiS^) is a positive 

definite matrix. 

Spectral Representation 

In the case of the damped harmonic oscillator, we saw explicitly that the dissipation 
was proportional to the coefficient of friction, 7. But for our quantum systems, the 
dynamics is entirely Hamiltonian: there is no friction. So what is giving rise to the 
dissipation? In fact, the answer to this can also be found in our analysis of the harmonic 
oscillator, for there we found that in the limit 7 — > 0, the dissipative part of the response 
function x" doesn't vanish but instead reduces to a pair of delta functions. Here we 
will show that a similar property holds for a general quantum system. 

We'll take the state of our quantum system to be described by a density matrix 
describing the canonical ensemble, p = e~@ H . Taking the Fourier transform of the 
Kubo formula (4.32) gives 

XiM = -ij^dt e-*Tr(e-^ 0,(0)]) 

We will need to use the fact that operators evolve as Oi(t) = C/ _1 Oi(0)C/ with U = e 
and will evaluate Xij( u ) by inserting complete basis of energy states 

POO 

Xij{uj) = ~i dt e iuJt J2 e ~ Em ^ [(m|O i |n)(n|e) j |m)e^- E ")* 

-{m\0 3 \n)(n\Oi\m)e i{En - Em)t ] 



-iHt 
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To ensure that the integral is convergent for t > 0, we replace uj — > uj + ie. Then 
performing the integral over J dt gives 

( O i ) mn ( ) nrrt ( Oj ) ran ( ^ i ) nm 



Xij(w + ie) = ^ 



e -EmP 



uj + E m - E n + ie uj + E n - E m + ie 

_ V^i)mn{Oj)nm ( e ~E m /3 _ e ~E„l3\ 

,,, _L R _ R _L ,V V ' / 



m,n 



which tells us that the response function has poles just below the real axis, 

uj = E n - E m - ie 

Of course, we knew on general grounds that the poles couldn't lie in the upper half- 
plane: we see that in a Hamiltonian system the poles lie essentially on the real axis (as 
e — > 0) at the values of the frequency that can excite the system from one energy level 
to another. In any finite quantum system, we have an isolated number of singularities. 

As in the case of the harmonic oscillator, in the limit e — > 0, the imaginary part of 
the response function doesn't disappear: instead it becomes a sum of delta function 
spikes 

m,n v ' m,n 

The expression above is appropriate for quantum systems with discrete energy levels. 
However, in infinite systems — and, in particular, in the quantum field theories that 
we turn to shortly — these spikes can merge into smooth functions and dissipative 
behaviour can occur for all values of the frequency. 

4.3.2 Fluctuation-Dissipation Theorem 

We have seen above that the imaginary part of the response function governs the 
dissipation in a system. Yet, the Kubo formula (4.32) tells us that the response formula 
can be written in terms of a two-point correlation function in the quantum theory. 
And we know that such two-point functions provide a measure of the variance, or 
fluctuations, in the system. This is the essence of the fluctuation-dissipation theorem 
which we'll now make more precise. 

First, the form of the correlation function in (4.32) — with the commutator and 
funny theta term — isn't the simplest kind of correlation we could image. The more 
basic correlation function is simply 

S^t) EE (0^)0,(0)) 
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where we have used time translational invariance to set the time at which Oj is evalu- 
ated to zero. The Fourier transform of this correlation function is 



Sij 



(u) = J dt e^Sijit) (4.34) 



The content of the fluctuation-dissipation theorem is to relate the dissipative part of 
the response function to the fluctuations S(cu) in the vacuum state which, at finite 
temperature, means the canonical ensemble p = e~$ H . 

There is a fairly pedestrian proof of the theorem using spectral decomposition (i.e. 
inserting a complete basis of energy eigenstates as we did in the previous section). Here 
we instead give a somewhat slicker proof although, as we will see, it requires us to do 
something fishy somewhere. We proceed by writing an expression for the dissipative 
part of the response function using the Kubo formula (4.32), 

= -\[xiAt)-xa(-t)\ 



2 



+\o(-t) (Oji-tWO)) - (Oi(p)O s (-t)) 



By time translational invariance, we know that (Oj(O)Oi(t)} = (Oj(—t)Oi(0)). This 
means that the step functions arrange themselves to give 6(t) + 0(—t) = 1, leaving 

4w = -\m)o 3 m + \(oA-t)om) (4.35) 

But we can re-order the operators in the last term. To do this, we need to be sitting 
in the canonical ensemble, so that the expectation value is computed with respect to 
the Boltzmann density matrix. We then have 

(Oji-^OiiO)) = Tre- pH Oj(-t)Oi(0) 

= Tre- fiH Oj(-t)e pH e- fiH Oi(0) 
= Tr e-^Oi^Oji-t + iP) 
= (0^-1(3)0^)) 

The third line above is where we've done something slippery: we've treated the density 
matrix p = e^^ H as a time evolution operator, but one which evolves the operator in 
the imaginary time direction! In the final line we've used time translational invariance, 
now both in real and imaginary time directions. While this may look dodgy, we can 
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turn it into something more palatable by taking the Fourier transform. The dissipative 
part of the response function can be written in terms of correlation functions as 



= -\ [(Oi(t)Oj(o)) - (Oi(t - i(5)OM) 



(4.36) 



Taking the Fourier transform then gives us our final expression: 



(4.37) 



This is the fluctuation-dissipation theorem, relating the fluctuations in frequency space, 
captured by S(u), to the dissipation, captured by x?'(u>). Indeed, a similar relationship 
holds already in classical physics; the most famous example is the Einstein relation 
that we met in Section 3.1.3. 

The physics behind (4.37) is highlighted a little better if we invert the equation. We 
can write 



where ub{oj) = (e^ — l)^ 1 is the Bose-Einstein distribution function. Here we see 
explicitly the two contributions to the fluctuations: the ub{oj) factor is due to thermal 
effects; the can be thought of as due to inherently quantum fluctuations. As usual, 
the classical limit occurs for high temperatures with f3u <C 1 where ns{oo) ~ k^T juj. 
In this regime, the fluctuation dissipation theorem reduces to its classical counterpart 



4.4 Response in Quantum Field Theory 

We end these lectures by describing how response theory can be used to compute some 
of the transport properties that we've encountered in previous sections. To do this, 
we work with quantum field theory 10 , where the operators become functions of space 
and time, 0(x, t). In the context of condensed matter, this is the right framework 
to describe many-body physics. In the context of particle physics, this is the right 
framework to describe everything. 

10 See http://www.damtp.cam.ac.uk/user/tong/qft.html for an introductory course on quantum field 
theory. 



Sij(u) 



2 [n B M + 1] 4-m 
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Suppose that you take a quantum field theory, place it in a state with a finite amount 
of stuff (whatever that stuff is) and heat it up. What is the right description of the 
resulting dynamics? From our earlier discussion, we know the answer: the low-energy 
excitations of the system are described by hydrodynamics, simply because this is the 
universal description that applies to everything. (Actually, we're brushing over some- 
thing here: the exact form of the hydrodynamics depends on the symmetries of the 
theory, both broken and unbroken). All that remains is to identify the transport co- 
efficients, such as viscosity and thermal conductivity, that arise in the hydrodynamic 
equations. But how to do that starting from the quantum field? 

The answer to this question lies in the machinery of linear response that we developed 
above. For a quantum field, we again add source terms to the action, now of the form 

#sourcc(t) = J ' d d ~ l X (f>i(S,t)Oi(S,t) (4.38) 

The response function x is again defined to be the change of the expectation values of 
O in the presence of the source 0, 

W , t ,)./^^^,0#,0 (4-39) 

All the properties of the response function that we derived previously also hold in the 
context of quantum field theory. Indeed, for the most part, the label x and rr'can be 
treated in the same way as the label Going through the steps leading to the Kubo 
formula (4.32), we now find 

XtJ (x,x>;t-t>) = - 0<[Oi(x, ^.^(f ,*')]) (4-40) 

If you've taken a first course on quantum field theory, then you know that the two-point 
functions are Green's functions. Usually, when thinking about scattering amplitudes, 
we work with time-ordered (Feynman) correlation functions that are relevant for build- 
ing perturbation theory. Here, we interested in the retarded correlation functions, 
characterised by the presence of the step function sitting in front of (4.40). 

Finally, if the system exhibits translational invariance in both space and time, then 
the response function depends only on the differences t — t' and x — x'. In this situation 
it is useful to work in momentum and frequency space, so that the (4.39) becomes 



— * — * — * 

6{Oi(k,u) = Xij(k,uj)(l) j (k,uj) 



(4.41) 



Electrical Conductivity 

Consider a quantum field theory with a U(l) global symmetry. By Noether's theorem, 
there is an associated conserved current J M = (J°, J 1 ), obeying <9 M J M = 0. This current 
is an example of a composite operator. It couples to a source which is a gauge field 



-^source 



J d d - x x (4.42) 



Here is the background gauge field of electromagnetism. However, for the purposes 
of our discussion, we do not take A^ to have dynamics of its own. Instead, we treat it 
as a fixed source, under our control. 

There is, however, a slight subtlety. In the presence of the background gauge field, 
the current itself may be altered so that it depends on A^. A simple, well known, 
example of this occurs for a free, relativistic, complex scalar field tp. The conserved 
current in the presence of the background field is given by 

= ie^dPy - (dVV] - e 2 AV</? (4-43) 

where e is the electric charge. With this definition, the Lagrangian can be written in 
terms of covariant derivatives V^p = d^ip — ieA^p, 



L = J d d ~ l x |<9 M v?| 2 + A^ = J d d ~ l x \V^\ 



(4.44) 



For non-relativistic fields (either bosons or fermions), similar A^ terms arise in the 
current for the spatial components. 

We want to derive the response of the system to a background electric field. Which, 
in more basic language, means that we want to derive Ohm's law in our quantum field 
theory. This is 

(Ji(k, u)) = aij(k, u)Ej(k, u) (4.45) 

Here is the background electric field in Fourier space and is the conductivity 
tensor. In a system with rotational and parity invariance (which, typically means in 
the absence of a magnetic field) we have cr^ = aSij, so that the current is parallel to 
the applied electric field. Here we will work with the more general case. Our goal is to 
get an expression for aij in terms of correlation functions in the field theory. Applying 
(4.41) with the perturbation (4.42), we have 

S(J,) = (J„) - (J M ) = -i I dt'd 3 x ([J^x,t),J u (x',t')})oA u (x',t') (4.46) 

J — oo 
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The subscript here means the quantum average in the state = before we turn 
on the background field. Let's start by looking at the term (Jj)o- You might think that 
there are no currents before we turn on the background field. But, in fact, the extra 
term in (4.43) gives a contribution even if - as we'll assume - the unperturbed state 
has no currents. This contribution is 

(Ji)o = e 2 A(V>V)o = eA lP 

where p is the background charge density. Notice it is not correct to set A, = in this 
expression; the subscript only means that we are evaluating the expectation value in 
the Ai = quantum state. 

Let's now deal with the right-hand side of (4.46). If we work in A = gauge (where 
things are simplest), the electric field is given by E-i = —Ai. In Fourier transform space, 
this becomes 

AM = ^1 (4.47) 

too 

We can now simply Fourier transform (4.46) to get it in the form of Ohm's law (4.45). 
The conductivity tensor has two contributions: the first from the background charge 
density; the second from the retarded Green's function 

„ a = + 2S&0 (448) 

IUJ IUO 

with the Fourier transform of the retarded Green's function given in terms of the 
current-current correlation function 

/oo 
dt d 3 x 6(t) e'O"*-*-*) ( [Jj (x , t) , Jj (0,0)]) 
-oo 

This is the Kubo formula for conductivity. 
Viscosity 

We already saw in Section 2 that viscosity is associated to the transport of momentum. 
And, just as for electric charge, momentum is conserved. For field theories that are 
invariant under space and time translations, Noether's theorem gives rise to four cur- 
rents, associated to energy and momentum conservation. These are usually packaged 
together into the stress-energy tensor T^, obeying d^T^ = 0. (We already met this 
object in a slightly different guise in Section 2, where the spatial components appeared 
as the pressure tensor and the temporal components as the overall velocity uA. 
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The computation of viscosity in the framework of quantum field theory is entirely 
analogous to the computation of electrical conductivity. The electric current is simply 
replaced by the momentum current. Indeed, as we already saw in Section 2.5.3, the 
viscosity tells us the ease with which momentum in, say, the x-direction can be trans- 
ported in the z-direction. For such a set-up, the relevant component of the current is 
T xz . The analog of the formula for electrical conductivity can be re- interpreted as a 
formula for viscosity. There are two differences. Firstly, there is no background charge 
density. Secondly, the viscosity is for a constant force, meaning that we should take 
the ui — > and k — > limit of our equation. We have 

/oo 
dt d 3 x 6(t) e^-^ ([T xz (x, t),T xz (6, 0)]> 
-oo 

and 

Xxz,xz(0,u) 

7] = lim 

iuj 

This is the Kubo formula for viscosity. 
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